Slippage of Newtonian liquids: influence on the dynamics of dewetting thin films 
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Slippage of Newtonian liquids in the presence of a solid substrate is a newly found phenomenon the 
origin of which is still under debate. In this paper, we present a new analysis method to extract the 
slip length. Enhancing the slip of liquids is an important issue for microfluidic devices that demand 
for high throughput at low pumping power. We study the velocity of short-chained liquid polystyrene 
(PS) films dewetting from non-wettable solid substrates. We show how the dynamics of dewetting is 
influenced by slippage and we compare the results of two types of substrates that give rise to different 
slip lengths. As substrates, Si wafers were used that have been coated by octadecyltrichlorosilane 
(OTS) or dodecyltrichlorosilane (DTS), respectively. Our results demonstrate that the dewetting 
velocity for PS films on DTS is significantly larger than on OTS and that this difference originates 
from the different slip lengths of the liquid on top of the two surfaces. For PS films of thicknesses 
between 130 nm and 230 nm we find slip lengths between 400 nm and 600 nm, depending on 
substrate and temperature. 



I. INTRODUCTION 

Thin films of liquid polymers play an important role in 
numerous technological processes ranging from lithogra- 
phy to biological membranes [lj . A key to understanding 
the stability of, e.g;, liquid coatings on solid surfaces was 
the calculation 0,111, Si IE 0] and the experimental deriva- 
tion of the effective interface potential [fl. M. l9l. [ToL [ill [l2| . 
In recent years, one major focus of interest has been 
to understand the dynamics and morphology of poly- 
mer films dewetting from non-wettable substrates [131 

11^ [IE, b^a IIS. IIE. 113 - The 

dewetting rate of the liquid hereby strongly depends on 
the flow velocity at the solid/liquid interface. Hitherto it 
was common sense that the occurrence of a non-zero flow 
velocity (slippage) at the interface is restricted to high 
molecular weight polymer films [28| . Only recently, new 
experimental techniques have revealed that also Newto- 
nian liquids may exhibit slippage [IE HE HH [HJ . The oc- 
currence and the nature of slippage is of large technolog- 
ical interest, since a sliding fluid can flow faster through, 
e.g., microfluidic devices. 

Many studies have focused on techniques to measure 
slippage. These approaches can be classified in direct 
and indirect methods to determine the fluid velocity pro- 
file near the solid/liquid interface. Direct methods use, 
e.g., tracer particles in combination with near-field laser 
velocimetry [IE HE] or fluorescence recovery after photo- 
bleaching techniques [IE HE H3| ■ Indirect methods aim 
at measuring the drainage force of an object moved in 
a liquid to calculate the amount of slippage. Common 
techniques use a surface forces apparatus [53, HH or an 
atomic force microscope (AFM) with a colloidal probe 
[3^ |. For details see, e.g., the review articles from Lauga 
et al. [13 or Neto et al. [H|. 

A common measure of slippage is the slip length, which 
can be understood as the length below the solid/liquid 
interface where the velocity extrapolates to zero. For 
simple fluids, the slip length is found to be independent 



of the shear rate [33j, [39( , but is influenced by the wetta- 
bility of the substrate as well as by the roughness of the 
wall: Molecular dynamics simulations for a simple liquid 
with 140° contact angle on the surface of interest found 
a slip length in the order of 30 diameters of the fluid 
molecules, but no-slip boundary condition for a wetting 
situation [40| . This result is in agreement with experi- 
mental data of hexadecane 0E [30, HH or glycerol [3l| . 
The role of surface roughness is ambiguous: simulations 
revealed that, depending on the pressure, surface rough- 
ness can increase or decrease the amount of slippage [41| ; 
experiments demonstrated the importance of the lateral 
length scale of the roughness [IE . 

In our studies, we investigate the dewetting process of 
thin liquid films, a typical scenario is depicted in Fig. [1] 
In order to extract the slip length, we develop a new 
method for the analysis of the dewetting rates. Dewet- 
ting can be described in three different stages [l6[ , start- 
ing from the very beginning of hole formation. In the 
first stage, the dewetted region grows exponentially in 
time and, in the vicinity of the dry area, a homogeneous 
thickening of the film can be observed; the hole does not 
yet exhibit a rim. In the second regime, the 'birth' of 
the rim takes place. In this stage, the rim has an asym- 
metric shape, and the radius grows linearly in time. In 
the third stage of dewetting, which starts at radii in the 
order of the slip length, surface tension rounds the rim 
which becomes more symmetric and grows from now on 
in a self-affine manner. Many experiments confirm this 
three-stage dewetting scenario [lj, [I3 [IE [2(| ■ 

In the third stage of a 'mature' rim, slippage is found 
to influence the dynamics in a substantial way 0, EH : 
For the case of a no-slip boundary condition or negligi- 
ble sliding friction at the solid/liquid interface, viscous 
dissipation (which occurs dominantly in the vicinity of 
the three phase contact line) results in a linear growth 
of holes; if, however, dissipation at the solid/liquid inter- 
face dominates the dewetting process, an R oc t 2 / 3 power 
law is expected. For the latter case of large slippage, the 
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third stage might be followed by a forth regime, when the 
rim height is large as compared to the slip length [l5[ . In 
this case, sliding dissipation again can be neglected and 
a linear growth R oc t is expected. 

In our studies we are interested in the stage of a ma- 
ture rim that grows in a self-affine manner in width and 
height. In the following we will demonstrate how we 
can deduce the slip length from the dynamics of the 
hole growth. We hereby restrict ourselves to liquids of 
Newtonian behavior and therefore investigate only the 
dewetting dynamics of short-chained polymer melts be- 
low the entanglement length. A short glance on Fig. [3] 
reveals that the growth law for the hole radius is clearly 
non-linear. We therefore suppose that slippage cannot 
be neglected. On the other hand, the polymer melt is 
not entangled. Hence, slippage is not expected to domi- 
nate the growth of holes at any stage. In the theoretical 
section we will describe a model assuming that neither 
viscous dissipation at the contact line nor sliding fric- 
tion at the solid/liquid interface entirely dominates the 
dynamics. We rather propose that the two different dis- 
sipation mechanisms are balanced and both have to be 
taken into account simultaneously [2l[ . 
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FIG. 1: Growth of a hole in a 130 nm thick PS(13.7k) film on 
DTS, captured by optical microscopy at 120 °C. 
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II. EXPERIMENTAL SECTION 

For the experiments, atactic polystyrene (PS, by Poly- 
mer Standard Service, Mainz, Germany) with a molec- 
ular weight of 13.7 kg/mol and a polydispersity of 
M w /M n = 1.03 was used. The samples were prepared 
by spin casting a toluene solution of PS onto mica, float- 
ing the films on Millipore™ water, and then transfer- 
ring them to hydrophobic substrates. For the hydropho- 
bization, we used two different self-assembled monolayer 
coatings on top of Si wafers (2.1 nm native oxide layer), 
octadecyltrichlorosilane (OTS) and the shorter dodecyl- 
trichlorosilane (DTS), respectively, prepared by standard 
techniques [iHf . By ellipsometry (EP 3 by Nanofilm, 
Gottingen, Germany), the thickness of the SAMs was 
found to be doTS = 2.4 nm and dors = 1-5 nm, re- 
spectively. Surface characterization by atomic force mi- 
croscopy (AFM, Multimode by Veeco, Santa Barbara, 
USA) revealed an RMS roughness of 0.09(1) nm (OTS) 
and 0.13(2) nm (DTS) at a (1 £tm) 2 scan size, and a static 
contact angle of polystyrene droplets of By = 67(3)° 
on both coatings. The surface tension of polystyrene is 
7 = 30.8 mN/m. 

The polystyrene films in this study are either 
130(5) nm or 230(5) nm thick. To induce dewetting, the 
films were heated above the glass transition temperature 
of PS (93 °C) to three different temperatures (110 °C, 
120 °C, 130 °C). After a few seconds, circular holes are 
born and instantly start to grow. An example of a typical 
temporal series captured by optical microscopy is shown 
in Fig. [T] Due to mass conservation of the liquid rims 
develop, surrounding each hole. Cross sections of such 
a growing rim are shown in Fig. [21 as measured by in 



FIG. 2: Snapshots of the rim profile of a dewetting PS(13.7k) 
film on OTS, measured in situ with tapping mode AFM. The 
dynamic contact angle stays constant during the growth of 
the hole. Here, the three phase contact line is shifted to the 
origin. 



situ AFM. The dynamic contact angle of the liquid is 
9dyn = 56(2)° and stays constant during dewetting. 

In Fig. [31 the optically measured radii of the emerging 
holes in the PS(13.7k) films as a function of time are 
depicted for two temperatures and the two types of silane 
substrate. Firstly, for each series, we find a non-linear 
growth law. To emphasize this, the data are additionally 
drawn in Fig. [¥] in a logarithmic diagram, revealing an 
algebraic growth of the radius in time with an exponent 
of a = 0.82(5). Secondly, dewetting progresses clearly 
faster on DTS than on OTS coated substrates. 

We qualitatively interpret this observation as follows. 
At the same temperature the polymer films on both sam- 
ples have identical properties: the viscosity r\ as well as 
the surface tension 7 do not depend on the substrate 
underneath. Additionally, the static contact angle 0y 
of polystyrene on both surfaces is constant within the 
experimental error. Therefore, the spreading coefficient 
S = -f(cos8 Y - 1) = -0.0188 N/m which is the driving 
force of the dewetting process is identical on OTS and 
on DTS substrates. The only parameter that could be 
different is the boundary condition at the solid/liquid in- 
terface of the two substrate. Since dewetting on DTS is 
much faster than on OTS, we expect a significantly larger 
slip length for PS(13.7k) on the DTS coating. 
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FIG. 3: Radius of dewetted region as a function of time, mea- 
sured for 130 nm thick PS(13.7k) films on OTS as well as on 
DTS substrates at two different temperatures. 
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FIG. 4: Data shown in Fig. [3] in logarithmic scale. For 
different dewetting temperatures and on the different sub- 
strates, the radii grow algebraically in time with an exponent 
of a = 0.82(5). 



III. THEORETICAL SECTION: ENERGY 
DISSIPATION MODELS 

To get more insight in the mechanisms involved in the 
dewetting process, a detailed analysis of the velocity is 
required. The velocity of the moving front V — R is given 
by the force balance between the driving force \S\ and the 
power P dissipated in the dewetting process, \S\ ■ V = P. 
Assuming no slip at the solid/liquid interface, energy is 
dissipated by viscous friction only within the liquid vol- 
ume and at the three phase contact line. Since the high- 
est shear rates arise in the vicinity of the contact line, 
viscous dissipation P v takes place dominantly at this line 
and is therefore independent of the size of the rim. How- 
ever, P v is influenced by the flow geometry near the con- 
tact line, i.e., by the dynamic contact angle 9d yn - Thus, 
the dewetting velocity V = v v is given by 



Vv — C v (0dyn) " — " 

V 



(1) 



with r] the viscosity of the liquid, and C v (9dyn) a function 
of the dynamic contact angle that counts for the flow 



geometry 

As shown in Fig. [2] the dynamic contact angle is tem- 
porally constant in our experiments. Thus, viscous dis- 
sipation does neither depend on hole size nor on time. 
On the other hand, the driving force 151 does not vary 
during the dewetting process. Hence, in the model of no- 
slip boundary condition, the dewetting velocity V = v v 
is constant, and the radius R grows linearly in time, 
R oct, which holds for Newtonian liquids as shown, e.g., 
in Ref. [H|. 

Since the data in Fig. [3] clearly do not show a linear 
growth of the radius i?, but rather a decreasing velocity, 
our experiments do not allow the assumption of a no-slip 
boundary condition. 

Another model, the so-called 'full-slip' or 'plug flow' 
model, introduced to explain the dewetting rate of en- 
tangled polymer melts, assumes large slip lengths at the 
solid/liquid interface [15| . Here, energy dissipation domi- 
nantly occurs at this interface over the distance of about 
the width of the rim, w. Hence, the power dissipated 
in the dewetting process is proportional to w, and force 
balance results in the velocity V = v Sl 
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with the slip length b. Self-similarity of the growing 
mature rim, (which means a constant dynamic contact 
angle, c.f. Fig. [2]), and mass conservation of the liquid 
yields w oc yR. Hence, in this model, the velocity de- 
creases with increasing hole size, v s oc I/VP, and the 
hole radius increases algebraically in time with an expo- 
nent a = 2/3, R oc i 2 / 3 . For simplicity, we can introduce 
the constant number K and write v s = Kj \[R. 

Although the data in Fig.|4]show an algebraic behavior 
R oc t a , the exponents fitted to the data achieve values 
clearly above 2/3; the mean value of the exponents of the 
four curves shown in Fig. [1] is a — 0.82(5). Hence, the 
model of full slippage also is not adequate to describe our 
experimental results. 

Another model combining energy dissipation at the 
three phase contact line and sliding dissipation at the 
solid/liquid interface is required. In literature there are 
some suggestions of a transition within the dewetting pro- 
cess, meaning that at an early stage viscous dissipation 
at the contact line dominates the dynamics of dewetting, 
R oc t, while in a later stag e the holes grow mostly due 
to slip effects, R oc t 2 / 3 [13 El El. In our experiments, 
for some small holes, we can see such a transition from 
an early stage of almost constant dewetting velocity to a 
later stage of decreasing velocity. Nevertheless, the model 
of slip-dominated dissipation is not appropriate for the 
experimental data, even not at the later stage. 

In our earlier study, Ref. [2l[, simple addition of 
both contributions of the dewetting velocity is suggested 
('combined model'), 
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FIG. 5: Except for very small holes, the data of the dewetting 
velocity V = R are in line when plotting versus 1/VR. Ex- 
trapolation for infinitely large holes yields to a finite velocity 

with v v = const, and v s = K/yR. Separation of the 
variables yields the function 

t-t = ^ ^i?-2 7 ^+ 2 7 2 ln^l + ^^ (4) 

7 = K/v v . 

Unfortunately, this implicit equation for the radius is 
rather cumbersome. Since it can be fitted to any data 
of hole growth, this function is not able to critically test 
the combined model and its assumptions. 



IV. TESTING THE COMBINED MODEL 

An alternative and more comprehensive way to check 
the combined model Eq. ([3]) is to plot the velocity versus 
1/^/R. If the combined model is valid, the data points 
should then appear on a straight line with intercept v v 
and slope K. As shown in Fig. [5j our experimental data 
indeed lie on such a straight line with a finite velocity 
Vq extrapolating 1/yR to zero, i.e., for infinitely large 
holes. This gives a first hint that the simple 'combined 
model' V = v v + v s is an adequate description for the 
dewetting rate of short chain polymer films. Just the 
data of small holes on DTS with radii below 4 /im deviate 
from this line. Since the rims of these holes are not yet 
in the mature regime, we can neglect these data points 
for further analysis. 

How can we further test that the combined model is 
indeed an appropriate description? 

If the intercept vq represents the dewetting velocity of 
viscous flow, v v , then it should fulfill three conditions: i) 
the temperature dependence of vq should be dominated 
by the inverse viscosity of the polymer melt, ii) vq should 
be independent of the initial film thickness /io, which is 
a very strong condition, and iii) the intercept vq should 
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FIG. 6: For evaluation of the combined model all measured 
data of dewetting velocities of PS(13.7k) films at 120 °C are 
plotted versus 1 / y/R- 
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FIG. 7: The finite velocity vo extrapolated for infinite large 
hole radii increases with increasing dewetting temperature. 
This behavior qualitatively compares well to the temperature 
dependence of the inverse melt viscosity. Additionally, vo 
does neither depend on the initial film thickness nor on the 
substrate coating. 



only depend on the contact angle, but not on the slip 
length. 

To test these three conditions, we repeated our dewet- 
ting experiments at three different temperatures and with 
films of two different thicknesses on both coatings DTS 
and OTS, cf. Fig. [5] for the experiments at 120 °C. As 
shown in Fig. the extrapolated velocity v increases 
for increasing dewetting temperature, and, moreover, the 
qualitative run compares well to the inverse viscosity of 
PS(13.7k) as measured by a rheometer. Note that vis- 
cous dissipation in the vicinity of the three phase con- 
tact line is not necessarily expected to exactly follow the 
1/n behavior, since there is a contribution additionally 
to the shear flow that derives directly from the contact 
line dynamics, i.e., thermally activated jumps over an 
energy barrier, pinning, etc. |44| . Thus, the first con- 
dition for vq concerning the temperature dependence is 
satisfied sufficiently. From Fig. [7] we also learn that the 
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initial film thickness ho does not systematically influence 
the extrapolated velocity v , although the measured ve- 
locity at finite hole size does depend on film thickness, 
cf. Fig. [HI With this, the second condition is fulfilled. 
Since the intercept vq is independent of initial film thick- 
ness, it also does not depend on the width of the rim 
and, consequently, not on the effect of slippage at the 
solid/liquid interface. Additionally, vq is not influenced 
by the substrate, OTS or DTS covered silicon wafers, on 
which the PS films are supposed to show different slip 
lengths. This satisfies the third condition, the indepen- 
dence of slippage of vq. Note that in general vo depends 
on the substrate, since viscous dissipation is determined 
by the dynamic contact angle Qdyn- In our experiments, 
however, we found the same 9d yn on both substrates and, 
consequently, same values for vq. With the described 
tests we found all three conditions satisfied and therefore 
can assume that the intercept vq indeed compares to v v , 
which allows further analysis of the data. 
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FIG. 8: Width of the rim as a function of the hole radius, ex- 
emplarily shown for a 130 nm thick PS(13.7k) film dewetting 
from OTS at 110 °C. By fitting Eq. © to the data we get 
C s = 1.96. 



V. RESULTS: THE SLIP LENGTH 

In a second step the velocity component v s is analyzed 
and the slip length b is extracted. In the case of full 
slippage at the solid/liquid interface the velocity is given 
by Eq. ([2]), where the width of the rim reads 



w = C s • ■ (5) 

Hence, the slope K shown in Fig. [5] is expected to be 

b 



3 V C s 



fha 



(0) 



First, the validity of Eq. is tested by plotting the 
width of the rim, w, versus the hole radius R, cf. Fig. [8] 
Here, w is taken as the lateral distance between the three 
phase contact line (x = 0) and the position where the rim 
height is dropped to 110% of the prepared film thick- 
ness ho, i.e., h(w) = 1.1 ho. Fitting Eq. |(SJ| to the 
data of 130 nm and 230 nm thick PS films, we obtain 
C s = 1.96(5) on OTS and C s = 2.1(1) on DTS covered 
wafers. 

Knowing C s , the spreading coefficient |S| = 
0.0188 N/m, the initial film thickness ho, and the vis- 
cosity t\ of the investigated films, cf. Fig. [JJ we can now 
use Eq. ([6]) to directly determine the slip length b from 
the slope K. Although the values for K are affected by 
the film thickness, cf. Figs. ® or [3Ji, this dependence is 
canceled out very well for the slip length, as shown in 
Fig. [9Jd. This again indicates that the combined model 
Eq. @ can be used to analyze the data of our experimen- 
tal system. We find b decreasing for increasing dewetting 
temperature. Additionally, the slip length of PS(13.7k) 
is clearly larger on DTS than on OTS, but the difference 
decreases with increasing temperature, a fact that cor- 
roborates the results for the slip lengths determined by 
rim shape analysis [45| . 
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FIG. 9: Left: Slope K of the straight lines in Fig. Has a 
function of the melt temperature. The parameter K depends 
on the initial film thickness as well as on the type of coating. 
Right: The slip length for polystyrene on top of both DTS 
and OTS covered Si wafers. The values are determined from 
K. 



VI. CONCLUSION 

We have developed a new analysis method to extract 
slip lengths from dewetting rates. The experiments give 
evidence that a model has to be put forward that com- 
bines frictional energy dissipation at the interface with 
viscous dissipation at the three phase contact line. We 
show that simple addition of two different velocity com- 
ponents is an adequate description for highly viscous 
Newtonian liquids dewetting from smooth surfaces. We 
moreover were able to demonstrate that this descrip- 
tion also captures a variable friction coefficient at the 
solid/liquid interface. The components v v and v s are 
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given by the well known expressions derived for the lim- 
iting cases no-slip and plug flow, respectively. From the 
slope of velocity data plotted versus 1 j\fR the slip length 
can be determined. 

For short-chained polystyrene melts on OTS and DTS 
covered wafers we found slip lengths between 400 nm and 
6 /zm. Slippage in this system decreases for increasing 
dewetting temperature. On the DTS coating the slip 
lengths are up to one order of magnitude larger than on 
OTS. 

So far, we have only hypotheses for the molecular 
mechanisms at the solid/liquid interface that give rise 
to the different slip lengths on OTS and DTS covered 
wafers. The commonly investigated system parameters 
like interaction forces between solid and liquid, i.e., con- 
tact angle and long-ranged dispersion forces, or substrate 



roughness are identical on both types of substrate. Thus, 
we have no straightforward explanation. The molecular 
origin of slippage in our system has to be postponed to 
future research. The simple analysis method of dewetting 
rates, however, is a powerful tool that allows for exten- 
sive studies of various systems to get a comprehensive 
picture of slippage of simple and complex liquids. 
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